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Influence of Homeotropic Anchoring Walls upon Nematic and Smectic Phases 

Masashi Torikai * and Mamoru Yamashita t 
Department of Physics Engineering, Faculty of Engineering, Mie University, 1515 Karnihama, Tsu 514-8507 

McMillan liquid crystal model sandwiched between strong homeotropic anchoring walls is 
studied. Phase transitions between isotropic, nematic, and smectic A phases are investigated 
for wide ranges of an interaction parameter and of the system thickness. It is confirmed that 
the anchoring walls induce an increase in transition temperatures, dissappearance of phase 
transitions, and an appearance of non-spontaneous nematic phase. The similarity between 
influence of anchoring walls and that of external fields is discussed. 

KEYWORDS: liquid crystal, McMillan model, homeotropic anchoring, nematic phase, smectic A phase, 
non-spontaneous nematic phase 
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1. Introduction 

Liquid crystals are very sensitive to conditions of 
boundary surfaces. The anchoring conditions at the 
boundary walls strongly influence the liquid crys- 
talline orderings if the systems are thin enough. The 
isotropic(I)-nematic(N) phase transition is easily affected 
by the anchoring walls since the discontinuity at the TN 
phase transition is small. It is known that a thin nemato- 
gen system sandwiched between homeotropic anchoring 
walls exhibits the I-N transition at higher temperature 
than the TN transition temperature of the system with- 
out boundary anchoring. 1-3 As the system becomes thin, 
the discontinuity at the TN transition decreases; eventu- 
ally the TN transition vanishes if the system is thinner 
than a critical thickness. 1 It is also shown theoretically 
that, even if the system exhibits I-smectic A(A) phase 
transition in bulk systems(i.e., in the absence of the 
walls), the homeotropic anchoring walls can induce a sur- 
face N phase; 4 such a N phase is called a non-spontaneous 
N phase. Corresponding to this result, strong N order in- 
duced near the anchoring walls has found in decylcyano- 
biphcnyl (10CB), which exhibits direct I-A transition in 
bulk systems. 5 The fact that such a wall-induced N or- 
der has not found in undecylcyanobiphenyl (11CB) and 
dodecylcyanobiphenyl (12CB) indicates that 10CB has 
stronger tendency to exhibit the N phase than others. 
Such a tendency is to be expected since the 10CB is 
the shortest molecule among the homologous series of 
cyanobiphenyl nCB exhibiting direct I-A transition in 
bulk systems. 

Although studies on a homologous series would provide 
a unified understanding of the phase behavior, there are 
no theoretical studies on the influence of anchoring walls 
upon members of homologous series. In the present pa- 
per, we investigate such influence in systematic way; in 
order to do so, we use a so-called McMillan model sand- 
wiched between strong homeotropic anchoring walls. The 
McMillan model can mimic several materials of a homol- 
ogous series by choosing a parameter a, contained in the 
model; depending on the a, the McMillan model exhibits 
three phases (I, N and A) or two phases (I and A). We re- 
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port the influence of the anchoring walls, i.e., the changes 
in phase transition temperatures, the critical thicknesses 
and their dependence on a, and the non-spontaneous N 
phase induced by anchoring walls. In the last section, we 
will discuss the similarity between the influence of an- 
choring walls and of the external fields on the McMillan 
model. 

2. Discrete McMillan Model 

Self-consistent equations for the discrete McMillan 
model in inhomogeneous systems have been derived 
in ref. 6, as a generalization of the original McMillan 
model. 7 Thus in this section we will shortly give defini- 
tions and formulae. 

The N and A phase order parameters are functions of 
the distance from one of walls in this system. We divide 
the system into layers parallel to the walls and assume 
the order parameters within each layer are constant. Let 
the number of layers be N + 2, and zeroth and (N + l)-th 
layers be in contact with walls. In the discrete McMillan 
model, the number of layers N can be interpreted as 
the system thickness. The N order parameter s n and A 
order parameter a n of n-th (0 < n < N + 1) layer are, 
by assuming that the z-axis is parallel to the director, 
defined as 
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where Pi(x) = (3x 2 — l)/2; 9 and z are the polar angle 
and z-coordinate of a molecule, respectively. The brack- 
ets (•)„ denote the thermal average over molecules in n-th 
layer. 

We take into account only the molecular fields pro- 
duced by neighboring molecules. Then the one-particle 
potential for a molecule in the n-th. layer is 



v n = -[Vs n + hW]P 2 (cos6) 

~ [aVa n + fi4 n) ] cos 



^z) />,(,-<* r>). (2) 



where Vs n and aVa n are molecular fields of the ro- 
th layer; hi^ = V'(s n _i + s n+ i — 2s„) and — 
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aV'(a n --i + <J n +i — 2<7„) arc molecular fields of neigh- 
boring layers, that vanish in homogeneous systems. The 
parameter a is an interaction strength for A phase; a can 
be assumed to be an increasing function of the molecu- 
lar length. 7 Using the discrete McMillan model potential 
(2), we obtain the self-consistent equations: 

h^),f3(aVa n + h^)), 



s n = I(fi{Vs n 
a n = J(/3(Vs r 



where functions I and J are defined as 
I( V ,() = -^lnZ(r 1 ,0, 

J( V ,0 = ^nZ(r 1 ,0, 
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Among the sets of solutions of eqs. (3), the thermo- 
dynamic stable set gives the minimum of a function 
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where the set {s n , <r n } is one of solutions to eqs. (3). The 
minimum of F is the thermodynamic free energy. 

3. Results 

For later discussions, it is suitable here to refer to the 
phase diagram of the McMillan model without bound- 
aries. 7 This model exhibits some different behaviors de- 
pending on the value of a. For a > 0.98, A phase directly 
melts into I phase; this transition is first-order. The N 
phase appears as an intermediate phase between I and 
A phases for a < 0.98. The FN transition temperature 
for a < 0.98 is T = 0.2202 (in units of V), which does 
not depend on a; the triple point of I, N and A phase is 
thus atripic = 0.98 and r tri pic = 0.2202. The FN phase 
transition is first-order for any a, while the N-A phase 
transition is first-order for 0.70 < a < 0.98 and second- 
order for a < 0.70; i.e., there is a tricritical point on the 
N-A coexisting line (at a tc = 0.70 and T tc = 0.1910). 

In the presence of walls or fields, the N order is induced 
even if the temperature is high; such a high temperature 
phase with finite N order is called a para-nematic phase. 
However, for simplicity, we do not use the term para- 
nematic but use the terms of systems without walls and 
fields; we call the high temperature phase I phase, the 
low temperature phase A phase, and the intermediate 
phase N phase, if it exists. 




0.222 



Fig. 1. Temperature dependence of order parameters of the n-th 
layer in a system of N = 18 and a = 1.0. From top to bottom, 
curves are for layers n = 2, 5, and 9. 



In this work, we assume the homeotropic anchoring 
condition, i.e., the director is perpendicular to the walls. 
In addition, we assume that the anchoring is so strong 
that the order parameters at the walls are extremely 
large. Thus we will solve the self-consistent equations 
under the boundary conditions s = s^+i = 1. Fur- 
thermore, we also assume that the strong A ordering is 
induced at the walls, i.e., Co = <Jn+i = F We set the 
potential parameters V = 1 and V — V/6. 

A typical behavior of order parameters of layers is 
shown in fig. 1 for a system of a = 1.0 and thickness 
N = 18. In our calculations, as shown in fig. 1, all the 
layers undergo each phase transition at the same tem- 
perature, i.e., the transition temperatures do not depend 
on the distance from the walls. Thus in order to find 
each transition temperature we can use averaged order 
parameters, (s) = J2 n S ™/ N and = J2 n a n/N, over 
the whole system. The transition at the same tempera- 
ture throughout the system is possibly a characteristic of 
systems under the complete wetting (by ordered phase) 
condition at the boundaries. In nematic liquid crystals, 
the phase transition of the whole system at the same tem- 
perature has observed in systems under complete wetting 
condition; 1 while under the incomplete wetting condition 
these systems undergo the "boundary-layer phase transi- 
tion" that occurs at a temperature higher than the bulk 
transition temperature. 8 

Figures 2 show the temperature dependence of or- 
der parameters (s) and (a) for a = 1.0. These figures 
clearly exhibit three important features: the increase in 
the transition temperatures, the existence of the non- 
spontaneous N phase, and the existence of critical thick- 
nesses below which the phase transitions disappear. The 
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Fig. 2. Temperature dependence of order parameters (s) and (a) 
for systems with a = 1.0. Curves correspond to the systems with 
thicknesses, reading from top to bottom, N = 7, 13, 18, 25, and 
infinity (i.e., bulk system). 



system with a = 1.0 exhibits direct I-A transition in the 
absence of the walls. By the influence of the walls the 
discontinuity at the transition decreases but there is no 
qualitative difference in systems thicker than N = 24. 
If the system becomes thinner than N = 24, the N 
phase appears between I and A phases. Thus in systems 
with thickness 14 < N < 24 three phases appear; i.e., 
the non-spontaneous N phase appears in these systems. 
As the system becomes thin, the discontinuities at the 
phase transitions decreases and I-N transition vanishes 
at N — 13; eventually the I-A transition also vanishes at 
N = 7. 

The critical thicknesses depend on the parameter a; 
the critical thicknesses for I-N and N-A phase transitions 
(ATj c N and N^ A , respectively) as functions of a are shown 
in fig. 3. The a-dependence of N^ A is strong while that 



of Nf N is weak, since the parameter a directly couples 
to the A phase order parameter in the one-particle po- 
tential (2). The N-A critical thickness seems to diverge 
as a approaches to 0.70 from above. This is because the 
discontinuity at the N-A transition is small near the tri- 
critical point, and thus the boundary anchoring condition 
strongly influences to the N-A transition. 

Figure 3 also shows that the curve Nf N terminates at 
a = 1.02, above which the intermediate N phase never 
appears no matter what the thickness is and only the I- 
A phase transition remains. Since the N phase does not 
appear in the bulk systems with a > atopic = 0.98, the 
non-spontaneous N phase appears in the systems with 




Fig. 3. The smectic A interaction parameter a dependence of 
critical thicknesses Nf N and N^ A . The circles and crosses de- 
note 7Vj c N and N^ A , respectively. Nf N (a) terminates at a = 1.02 
(indicated by a solid circle), since the nematic phase disappears 
if a > 1.02. 



^triple < a < 1.02. 

Such a non-spontaneous N phase just above the atripie 
corresponds to experimental results as follows. The a vs 
temperature phase diagram for the McMillan model in 
the absence of boundaries and fields can be compared 
with the alkyl chain length vs temperature phase dia- 
gram because of the relation between a and the molecu- 
lar length. 7 These phase diagrams are very similar as 
shown in figs. 2 and 6 of ref. 7. The experimental re- 
sult for 10CB liquid crystals sandwiched between walls 5 
shows that the non-spontaneous N phase appears at the 
surface. Since the 10CB is the shortest molecule among 
members of the homologous series of nCB exhibiting I-A 
transition, the 10CB can be assumed to correspond to 
the material with the parameter a just above the triple 
point atripie- 

The thickness dependences of the transition temper- 
atures are shown in figs. 4(a) and (b) for a = 1.0 and 
a = 0.88, respectively. From these figures, it seems that 
the increments of the transition temperatures, AT(iV) = 
T(N) - T bulk , depend on the system thickness as AT ~ 
N~ 2 in a very thin region. Such dependences differ from 
the prediction of the Kelvin equation: 9 

2 7h 



AT(N) 



71 



^bulk 



(7) 



N L 

where L is the latent heat per volume; 7h and 71 are, re- 
spectively, the surface tensions (per a wall) of higher and 
lower temperature phases. According to the Kelvin equa- 
tion, the thickness dependence of AT(N) is AT(N) ~ 
iV _1 . This difference from the Kelvin equation clearly 
shows that the assumption on which the Kelvin equa- 
tion based is not valid in very thin systems. The Kelvin 
equation is derived by assuming that the latent heat and 
surface tensions are independent of the system thick- 
ness, and this assumption is accurate for sufficiently 
thick systems. In fact as shown in fig. 5, the entropy 
of transition (i.e., L/T) strongly depends on the thick- 
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Fig. 4. Phase diagrams for (a.) a = 1.0 and (b)a = 0.88 systems. 

(a) T bulk denotes the transition temperature without walls; AT 
is the change in transition temperature induced by anchoring 
walls. The horizontal axis is normalized by the critical thickness 
N c = 8. The letters I, N, and A denote isotropic, nematic, and 
smcctic A phases, respectively. Circles are the critical points. 

(b) The higher and lower curves denote I-N and N-A coexisting 
lines. The horizontal axis is normalized by the critical thickness 
N c = 13. Crosses denote the transition temperatures of systems 
without walls. 



ness in thin region; the Kelvin equation is valid for much 
thicker systems than the system we considered here. The 
Clapeyron-Clausius equation, which is valid irrespective 
of the system thickness, 2 is appropriate to treat such thin 
systems. 

4. Summary and Discussions 

In summary, we have investigated the influence of 
the strongly homeotropic anchoring walls on the McMil- 
lan model system. We have shown that the existence 
of the critical thickness for I-N, N-A and I-A transi- 
tions. This result is analogous to the existence of criti- 
cal thickness for I-N transition in Landau-de Gennes ne- 
matogen model. 1 We have found that the system with a 
just above the atopic exhibits non-spontaneous N phase 
induced by anchoring walls; such a behavior is consis- 
tent with the behavior of a homologous series nCB. 5 
We have also obtained the I-N, N-A and I-A transition 
temperatures as functions of the system thickness. We 
have confirmed that the Kelvin equation is not valid for 
very thin system; the transition temperature behaves as 
T(N) - T bulk - N~ 2 instead of the prediction of the 
Kelvin equation T(N) - T bulk - N' 1 . 




Fig. 5. Changes in entropy per volume at phase transitions of 
a = 1.0 system. The horizontal dashed line at AS/V = 1.72 
denotes the entropy change in the bulk system. The curves with 
o, x, and o correspond to AS/V of I-N, N-A, and I-A phase 
transitions, respectively. 



We emphasize that the effects of anchoring walls are 
qualitatively similar to the effects of the external field 
as indicated in the early studies. 1,2 Let us introduce two 
external fields h s and h a , conjugate to the order param- 
eters s and a, respectively. We note that the squares of 
electric and magnetic fields correspond to h s , while h a 
has no experimental counterpart; thus usually only the 
field h s has been considered in literatures, 10-17 and h a 
has been omitted. In the following, we first consider only 
the influence of the field h s . The increase of the transition 
temperature and the non-spontaneous N phase occur un- 
der the influence of h s ; in this sense, the influence of the 
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Fig. 6. The transition temperatures T(h s ,h,j) for a = 1.0 sys- 
tem as a function of external fields h s and h a . This figure shows 
field-dependence of AT/T , where AT = T(h s , K) -T(0, 0) and 
To = T(0, 0) is the transition temperature in the absence of ex- 
ternal fields; this figure is restricted to a plane h = h 3 = h a for 
simplicity. The curve in0</t//t c <0.1 denotes I-A transition 
temperature; higher branch and lower branch in h/h c > 0.1 de- 
note I-N and N-A transition temperatures, respectively. Circles 
indicate the critical points of the corresponding phase transi- 
tions. 
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field h s and of the anchoring walls are qualitatively sim- 
ilar. However, the similarity is incomplete with respect 
to the existence of the critical points. By the application 
of sufficiently strong field h s the I-N transition vanishes 
but the N-A transition does not, i.e., the critical point 
on the N-A coexisting line of h s -T phase diagram does 
not exist. The N-A transition in a bulk McMillan sys- 
tem changes from first-order to second-order transition 
by the influence of h s ; i.e. the tricritical point appears on 
the N-A coexisting line instead of the critical point. 10-12 
Thus we can conclude that the similarity is broken since 
both the I-N and N-A coexisting line terminate at critical 
points in the system under the influence of the anchor- 
ing walls, as shown in fig. 4. However, if we consider the 
external field h„ in addition to h s , the tricritical point 
on the N-A coexisting line changes to a critical point in 
the bulk system; 6 then the similarity is again recovered. 
In order to show this similarity explicitly, we calculate 
the increment in transition temperatures under the ap- 
plication of both the h s and h a (see ref. 6 for details of 
the transition temperatures under the external fields). 
We plot on fig. 6 the increment of the transition tem- 
perature AT = T{h Sl h a ) — T(0, 0), under a condition 
h s = h a for simplicity. The similarity between the influ- 



ence of anchoring walls and of external fields is clearly 
seen in fig. 4 and fig. 6. This fact shows that, in consid- 
ering the anchoring effects in the relationship with the 
external field, it is necessary to introduce the fictitious 
field h a together with h s . 
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